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The test masses in gravitational-wave detectors will be sensitive not only to astrophysical gravi- 
tational waves, but also to the fluctuating Newtonian gravitational forces of moving masses in the 
ground and air around the detector. These effects are often referred to as gravity gradient noise. 
This paper considers the effects of gravity gradients from density perturbations in the atmosphere, 
and from massive airborne objects near the detector. These have been discussed previously by 
Saulson, who considered the effects of background acoustic pressure waves and of massive objects 
moving smoothly past the interferometer; the gravity gradients he predicted would be too small to 
be of serious concern even for advanced interferometric gravitational-wave detectors. In this paper 
I revisit these phenomena, considering transient atmospheric shocks, and estimating the effects of 
sound waves or objects colliding with the ground or buildings around the test masses. I also consider 
another source of atmospheric density ffuctuations: temperature perturbations that are advected 
past the detector by the wind. I find that background acoustic noise and temperature fluctuations 
still produce gravity gradient noise that is below the noise floor even of advanced interferometric 
detectors, although temperature perturbations carried along non-laminar streamlines could produce 
noise that is within an order of magnitude of the projected noise floor at 10 Hz. A deflnitive 
study of this effect may require better models of the wind flow past a given instrument. I also find 
that transient Shockwaves in the atmosphere could potentially produce large spurious signals, with 
signal-to-noise ratios in the hundreds in an advanced interferometric detector. These signals could 
be vetoed by means of acoustic sensors outside of the buildings. Massive wind-borne objects such as 
tumbleweeds could also produce gravity gradient signals with signal-to-noise ratios in the hundreds 
if they collide with the interferometer buildings, so it may be necessary to build fences preventing 
such objects from approaching within about 30m of the test masses. 

PACS numbers: 04.80.Nn, 95.55.Ym 



I. INTRODUCTION 

Interferometric detectors such as LIGO and VIRGO 
rely on exquisite sensitivity to the positions of hanging 
test masses in order to detect the perturbations of passing 
gravitational radiation. The sensitivity is so great that 
the measurements can also be affected by fluctuations in 
the local Newtonian gravitational field, which create tiny 
accelerations of the mass. This noise source, known as 
gravity gradient noise or Newtonian gravitational noise, 
is caused by the near-field gravity of masses moving near 
the interferometer, and is not to be confused with the 
far-field propagating gravitational waves that the instru- 
ments are intended to measure. 

Gravity gradient noise has the potential to be quite 
insidious, since it cannot be shielded by improvements 
to the test-mass vibrational isolation. The only effec- 
tive way to eliminate gravity gradient noise is to elimi- 
nate the moving masses that create the perturbing fields. 
Fortunately, though, the strongest perturbations to the 
local gravitational field are at frequencies well below the 
detectors' pass-bands. Since the proposed terrestrial in- 
terferometric detectors all have sensitivity cutoffs around 
3 Hz or higher, we need only worry about noise sources 
that can perturb the local gravity field on timescales less 
than about 0.3 seconds. Most of the noise sources that 
I consider are motivated by the expected sensitivity of 
advanced LIGO interferometers, which were originally 
projected to have a low-frequency cutoff around 10 Hz, 



and instrumental noise oi Sh ^ 2 X lO-^^Hz"^ in the 
gravitational- wave signal output at that frequency ||l| . 
Although changes in instrumentation technology will 
modify the ultimate sensitivity goals of LIGO, this "stan- 
dard" advanced LIGO noise level is a good reference 
point when considering new noise sources. Also, as 
pointed out below, gravity gradient noise will make it 
difhcult to push the detector noise much below this level 
at 10 Hz, regardless of improvements to the interferome- 
ters. 

Saulson |^ was the first to estimate the effect of 
gravity gradient noise on terrestrial interferometric de- 
tectors, considering the effects of seismic waves pass- 
ing through the earth and of sound waves in the air. 
In both cases he found the spectral density of noise 
in the interferometer path-length difference to be less 
than 10~^^m^Hz~^ around 10 Hz, corresponding to 
noise in the gravitational-wave signal at levels less than 
10~'^^Hz~^ for a 4 km interferometer. By compari- 
son, this is significantly less than the noise floor of 
'--^ 2 X 10~^^Hz^^ that advanced LIGO interferometers 
expect to achieve at 10 Hz. More recently, a detailed 
analysis has been made of seismic gravity gradient 
noise; this study indicated that seismic gravity gradi- 
ent noise would be within a factor of 2 of the advanced 
LIGO noise floor at 10 Hz for most times, and could actu- 
ally exceed this noise floor during seismically noisy times, 
making seismic gravity gradients a significant barrier to 
improvements in low-frequency sensitivity. It therefore 
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seems prudent to revisit the issue of atmospheric gravity 
gradients as well. 

In Sec. II I consider gravity gradients caused by at- 
mospheric pressure perturbations — the same noise source 
considered by Saulson. Attention is paid, however, to 
the effects of the ground, and of buildings that reduce 
the pressure noise in the immediate vicinity of the in- 
terferometer test masses. Nonetheless, I find that these 
tend only to weaken the gravity gradient noise in the 
pass-bands of interferometric detectors, reinforcing the 
conclusion that this noise source is not of any great con- 
cern. 

A much stronger source of high-frequency density per- 
turbations in the atmosphere is the presence of tempera- 
ture fluctuations, which are advected past a detector by 
the wind. In Sec. Ill I analyze this as a potential source 
of gravity gradient noise. However, I find that while 
small-scale temperature perturbations can produce high- 
frequency temperature fluctuations at any given point, 
they do not produce the same high-frequency fluctua- 
tions in the test mass position, since a given pocket of 
warm or cool air will affect the test mass gravitation- 
ally over the entire time that it is in the vicinity of the 
test mass, which is on the order of seconds. This pro- 
duces a cutoff in the noise spectrum above a few tenths 
of a Hz. The presence of turbulent vortices around the 
interferometer buildings can increase the high-frequency 
component, but still probably not enough to show up in 
the gravitational-wave noise spectrum. 

In Sees. IV and ^ I turn away from sources of back- 
ground noise to consider possible sources of transient 
gravity gradient signals that might be detected as spuri- 
ous events in the gravitational-wave instruments. Sec. IV 



extends the analysis in Sec. II to look at the effects of at- 
mospheric Shockwaves, such as might be generated by an 
explosion or supersonic aircraft. I find that sources such 
as these can indeed produce detectable signals that might 
be interpreted spuriously as gravitational- wave events. 
However, such signals would easily be vetoed using acous- 
tic monitors outside the interferometer buildings. 

Sec. analyzes the gravity gradients produced by indi- 
vidual objects, such as wind-borne debris, moving around 
outside the interferometer buildings. Saulson considered 
this effect for the case of objects moving with fairly uni- 
form velocity, but typically, in order to produce signifi- 
cant signal above 3 Hz, an object's motion must be chang- 
ing on timescales of less than 0.3 seconds. In particular, I 
find that objects colliding with the interferometer build- 
ings produce much stronger signals than objects simply 
passing by the buildings. As an example, tumbleweeds 
at the Hanford LIGO facility will be a steady source of 
spurious signals in advanced detectors if they are allowed 
to collide with the end stations. Preventing such signals 
requires shielding a region of a few tens of meters around 
the end station, screening any wind-borne debris that 
masses more than a few hundred grammes. 



groups and possible directions for further research. 



II. ATMOSPHERIC PRESSURE WAVES 

Pressure perturbations are the only type of atmo- 
spheric gravity gradient noise considered by Saulson 
The derivation in this section gives largely the same re- 
sult as his. 

Consider a plane pressure wave with frequency / prop- 
agating through a homogeneous airspace at some sound 
speed c. If the fractional pressure change Sp/p is small, 
it will induce an adiabatic density change Sp/p — 6p/jp, 
where 7 « 1.4 is the ratio of heat capacities at constant 
pressure and constant temperature for air at normal tem- 
peratures. The gravitational acceleration produced by 
this wave in the direction of propagation is: 



GzSp 



dV 



2Gpc 
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Spit + 1/4/) , (1) 



Sec. V] presents some concluding remarks, including 
recommendations to the gravitational-wave experimental 



where p ~ 1.3kg m~^ and p « lO^Nm"^ are the ambient 
air density and pressure, and Sp{t) is the pressure pertur- 
bation measured at the same point as the acceleration is 
being measured. By symmetry, there is no acceleration 
transverse to the wave. 

Now consider sound waves in the vicinity of the inter- 
ferometer. First, since the interferometer is only sensitive 
to motions of the test mass parallel to the arms, the grav- 
itational acceleration is reduced by a factor cos 9, where 
9 is the angle between the propagation direction and the 
interferometer arm. 

Second, the interferometer test mass is inside a build- 
ing, which can in principle be used to suppress noise 
within a distance r,„in of the test mass. Roughly 
speaking, this results in a high-frequency cutoff factor 
C(27r/rinin/c), where the function C{x) depends on the 
precise shape of the building, the manner in which it 
refiects sound waves, and many other factors, but is nor- 
mally close to 1 for a; < 1. For instance, if one simply 
removes from the volume integral in Eq. ([^) a cylinder 
with length and diameter both 2r,„in aligned with the 
z-axis, then C{x) ~ 1 for a; < 1, but oscillates with an 
amplitude of ^ 0.3 for x ^ 1. This function is shown 
in Fig. |l|. The constant-amplitude oscillations of C'{x) 
for large x refiect the assumption that the sound wave 
has a coherence length much longer than the building 
size, so the field between the two ends of the excluded 
cylinder is generated almost entirely by the first half- 
wavelength beyond each cap. Realistically, the actual 
high-frequency behaviour of C{x) will depend on how 
the sound waves bend and scatter around the building; 
however, this should not change the order of magnitude 
of C{x). The behaviour shown in Fig. |l| is therefore prob- 
ably a good estimate of the true cutoff function for a; ^ 1, 
and a reasonable order of magnitude estimate for a; ^ 1. 
For the LIGO end stations, r^^i^ is of order 5 metres, 
giving X ^ //(lOHz); the factor will not be too far off 
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FIG. 1: Plot of the factor C{x) indicating the reduction in 
sonic gravity gradient noise due to setting the density pertur- 
bations to zero within a cylinder of diameter and length 2rmin 
centered on the test mass (and aligned with the wave), ver- 
sus X = 27rrmin/A where A is the wavelength. The oscillatory 
behaviour above x = 1 results from the ends of the cylinder 
coming in and out of phase with each other. A more accurate 
model of the scattering of sound waves off an interferometer 
end station building would modify the behaviour of C{x) for 
X > 1, but probably not change its magnitude significantly. 



for the frequencies of greatest interest. More precise esti- 
mates would depend on the specific architectural details 
of a particular facility. 

Third, the interferometer is on the ground, not in ho- 
mogeneous empty space. For simplicity I assume that 
the waves are almost entirely reflect ed of f the ground; 

In this case 



I will justify this assumption in Sec. 11 A 



the gravity gradient in directions parallel to the ground 
contributed by the reflected wavefront is the same as if 
the wavefront were extended below ground, while the 
pressure perturbations measured by detectors near the 
ground (much less than a wavelength) will be doubled. 
The acceleration experienced by an interferometer test 
mass is therefore: 



= ^cos(^)C(2^/r„ 



Jc)Sp{t+l/4f). (2) 



The gravitational wave signal h{t) in the interferom- 
eter is related to the acceleration of one of the test 
masses by h{t) — g{t)/L, or in frequency space h{f) — 
(27r/)~^g(/)/L, where L is the length of the interferom- 
eter arm. Thus: 
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Assuming stationary noise, the one-sided spectral density 



Sh{\f\) is given by {KfMf'Y) = SHi\f\)S{f- f), where 
(. . .) denotes an expectation over all random phases of 
all plane wave modes contributing to the noise, and * 
denotes complex conjugation. Taking mode amplitudes 
and directions to be uncorrelated, this gives: 
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where Sp{\f\) is the acoustic noise spectral density mea- 
sured outside the building in the vicinity of a particular 
test mass. Since the two test masses in an arm are many 
wavelengths apart, their noise will be uncorrelated, and 
will thus add in noise power. Similarly, the noise from 
the two arms will add in power. (Actually this is a bit 
of an overestimate, since the noise in the motion of the 
test masses at the corner station will be somewhat cor- 
related.) Noting that {cos^ 9) = 1/3, one finds that the 
total noise in the gravitational wave signal is: 
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where i denotes a particular test mass in the interfer- 
ometer, r^ln is the dead air radius about the ith test 
mass, and >S'p''*(|/|) is the acoustic noise spectrum mea- 
sured outside the building enclosing that test mass. 

Infrasound noise spectra should be taken at the actual 
interferometer sites, but one can make estimates based 
on typical terrestrial atmospheric noise. An empirical 
study Q collected 256 power spectra of 1-16 Hz infra- 
sound data from a rural forest 50km from New York City 
over a period of months, and found that the average noise 
spectrum Sp{f) was relatively flat at 6-16 Hz, though 
with widely varying amplitude: 25% of the spectra had 
noise under - lOOnbarVHz, 50% under - 300nbarVHz, 
and 75% under ^ lOOOnbar'^/Hz. I use Eq. (||) to com- 
pute the corresponding noise in the LIGO detector. As- 
suming that the end masses (with rmin 5m) dominate 
the contribution to the gravity gradient noise and con- 
tribute equally, the noise in the gravitational-wave signal 
around 10 Hz is: 



(6 X lO-^'^Hz" 
xC(//10Hz) . 
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(6) 



The results are plotted in Fig. ^, using infrasound power 
spectra read off of Fig. 3 of ^ . 

Even the third-quartile power spectrum is between 
two and three orders of magnitude below the expected 
noise floor of 2 x 10~''^Hz~^ at 10 Hz projected for ad- 
vanced LIGO interferometers. Thus ambient infrasound 
is probably a negligible effect for determining the noise 
floor for most interferometric gravitational-wave detec- 
tors. Nonetheless the issue cannot be completely resolved 
without infrasonic noise data from the actual interferom- 
eter sites. 
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FIG. 2; Plot of dimensionless strain noise \/fSh(J) ver- 
sus frequency / for infrasonic atmospheric gravity gradients. 
The solid curve is the projected noise floor for advanced LIGO 
detectors; the dotted curves are the first, second, and third 
quartiles of noise produced by gravity gradients from ambi- 
ent pressure waves. Data for the infrasonic noise power are 
taken from Fig. 3 of [Q. Clearly the pressure waves would 
not contribute significantly to LIGO noise in any frequency 
range. 



A. Ground absorption 

In the derivation above I treated the sound waves 
as being reflected off the ground. Now consider what 
happens if a sound wave is absorbed by the ground. 
The energy flux in a traveling compression wave is 
C^/^ p~^/^ {{Sx / x)"^) , where C is the compression mod- 
ulus of the medium (Cair = 7Pair), and {{Sx/xY) is the 
average squared dimensionless compression factor over a 
wave cycle. The gravity gradient noise induced by such a 
wave goes as c^(((5p)^), where {{Sp)^) ^ p{{Sx/x)^), and 
c — \J C I p is the wave speed in that medium. Therefore, 
if a sound wave is completely absorbed by the ground, 
the resulting ground motions will produce gravity gradi- 
ent noise contributions in the ratio: 
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Now for Cair = 7Pa 

lO^Nm-^, Pair = 1.3kg m-^^ Pground ~ 1.8 X lO^kgm- , 
the ratio turns out to be of order 10^; that is, if sound 
waves were completely absorbed by the ground, the re- 
sulting ground vibrations would produce gravity gradi- 
ent noise levels about 1000 times greater than the atmo- 
spheric gravity gradients. However, it was shown in |^ 
that seismic gravity gradients are only just large enough 



to worry about. So the only way that atmospheric grav- 
ity gradients can be larger or of the same order as seismic 
gravity gradients is if the waves are mostly reflected off 
of the ground. This was one of the assumptions used in 
deriving Eq. (||). 



III. ATMOSPHERIC TEMPERATURE 
PERTURBATIONS 



The largest small-scale atmospheric density perturba- 
tions are caused not by pressure waves but by tempera- 
ture perturbations. As heat is transported up through 
a convective atmospheric layer, convective turbulence 
mixes pockets of warm and cool air to form temperature 
perturbations on all lengthscales down to a few millime- 
tres. On the timescales of interest (less than a second) 
these perturbations are effectively "frozen" into the air- 
mass, while pressure differences disperse rapidly in the 
form of sound waves. Perturbations in the air density 
p (X p/T are therefore caused predominantly by the tem- 
perature perturbations, which are typically several orders 
of magnitude larger than the pressure perturbations. Al- 
though they are frozen into the airmass, these tempera- 
ture perturbations can cause rapid time-varying density 
fluctuations 5p — —p5T/T as the wind carries them past 
a point in space. This is the primary source of "seeing" 
noise that affects optical astronomy. 

The appendix to this chapter gives a rigorous mathe- 
matical derivation of the gravity gradient noise spectrum 
due to these temperature perturbations. This section 
gives a qualitative derivation that reproduces the final 
result to order of magnitude. 

The gravity gradient signal at some frequency / is 
caused by pockets of warm or cool air with some length- 
scale I being advected past the interferometer test mass 
at a speed w, where I ~ vjlix] . Consider a single such 
pocket of air with a temperature perturbation bT away 
from the ambient temperature T. The gravitational ac- 
celeration produced in the instrument as a function of 
time t is g^(t) = Gpl^{6T/T)x{t)r-^{t), where p is the 
ambient air density, r{t) is the distance of the air pocket 
from the test mass as it is blown past, and x{t) is this dis- 
tance projected onto the axis of the interferometer arm. 
This geometry is sketched in Fig. ||. 

Now in general, the noise power spectral density in any 
quantity a due to a background of independent, uncor- 
related events is 5'a(|/|) = (2/At)|a(/)p, where d{f) is 
the Fourier spectrum from a single event and At is the 
spacing between events. Assuming uncorrelated pockets 
of air, then independent pockets of air arrive along any 
given streamline at intervals At ~ Z/v, and streamlines 
separated by more than I add noise incoherently (i.e., 
add linearly in power). This gives the following noise 
spectrum: 
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FIG. 3: Schematic of a pocket of air with temperature perturbation AT over a lengthscale /, being advected past a test mass 
k along a streamhne 5'. The pocket has an instantaneous velocity v(t) and position r{t) relative to the test mass, and x{t) is 
the projection of r{t) onto the axis of the interferometer arm. {S} denotes a reference plane intersecting orthogonally with all 
streamlines S. 



where /jsj denotes an integral over a plane crossing all 

streamlines 5', and Gs{f) is the Fourier transform of the 
function G{t) = x{t)/r^{t) taken along a given stream- 
line. The quantity ST'^{1) is the average squared tem- 
perature difference between points a distance I apart. 
Turbulent mixing theory, as well as actual micromete- 
orological measurements, predict a power-law behaviour 
for small separations: ST'^{1) ^ , where p is typi- 
cally 2/3. This applies for horizontal separations I up to 
of order 50 times the height of a given air pocket above 
the ground For streamlines more than a metre or so 
above ground level, then, this behaviour for 5T^ should 
be good out to distances I ^ 50m, corresponding to fre- 
quencies ^ 0.2Hz(i;/10ms^^). 

Eq. gives the gravity gradient noise on a given 
test mass in the interferometer. Denoting the test 
masses by the index k = 1 ... 4, and assuming that 
each test mass contributes independently to the noise 
in the gravitational- wave signal h, one has S';i(|/|) = 
(27r/)-4L-2^^S'g(|/|). Combining this with Eq. §) 
and the relation I ~ v/iirf yields: 



xj]/ dA\Gs.Af)\' 



(9) 



The more rigorous analysis in the appendix (Sec. ^) 
gives an expression with roughly the same form, but cov- 
ers the factors of order unity, and also accounts for the 
fact that wind speed can vary along a streamline and 



between streamlines. The more accurate formula is 

2 

LT 



Sum ^ 2^[%] 4(27r/)-(P+7)sin(p7r/2)r(p + 2) 



y-Y. Fs.kifrGsAf)wdA, (10) 

where w is the wind speed of the streamline as it crosses 
the plane of integration, and Fs.kif) and Gs.kif) are 
Fourier transforms of functions Fs.k{t) and Gs,k{t) de- 
scribing the motion of a point along a streamline S past 
a test mass fc, of the form: 
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(11) 
(12) 



It is worth noting that the frequency structure of 
Sh{\f\) depends on the time behaviour of the functions 
x{t) and r(t) describing the position of a point on a 
streamline relative to a test mass, which can be some 
distance away. The minimum distance rmin from the test 
mass to the passing air is thus an additional important 
scale in the problem: if x{t) and r{t) change significantly 
only on timescales ~ rmin/f, then the noise spectrum will 
be cut off at frequencies ^ v/2i:r-axin- 

By comparison, the temperature noise spectrum 
S't(|/|) measured at a point depends only on the lo- 
cal properties of the atmosphere at that point. Ap- 
plying similar order-of-magnitude arguments, one can 
write S't(|/|) {2l/v)cllP\H{f)\^, where H[t) - 1 
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for \t\ < l/v and otherwise. At high frequencies 
> 0.2Hz(t!/10ms~^) the system involves only one length- 
scale I, giving the spectrum a power-law dependence: 

Srilfl) ~ 24«^'(27r/)'(f+i) . (13) 

This is the same, to order of magnitude, as the exact 
result given in Eq. (A3). 



A. Uniform airflow 

The gravity gradient noise is easy to compute from 
Eq. (|l^ for the case of uniform airflow parallel to the 
ground with some constant velocity v. Placing the ref- 
erence plane orthogonal to v and passing through the 
test mass, the equations of motion for a given streamline 
past the test mass become quite simple: v{t) ^ w ~ v, 
x{t) = vtcos6'-|-ro sin^cos V', r{t) = (rp -l-w^t^)^/-^, where 
To is the distance from the test mass to the nearest point 
on the streamline, and rg cos "0 is the projection of this 
distance onto the ground. The geometry of the airflow is 
shown in Fig. ^. It is easy to show that: 

G{f) = -^[tcos0Ko{27rfro/v) 

+ sm9cosij Kii2TTfra/v)] , (14) 

Fif) = vP+'G{f) , (15) 

where Kq and Ki are moddified Bessel functions of the 
second kind of orders and 1. I perform the integral over 
the above-ground half of the reference plane, out from 
some radius that is roughly the closest distance that 
the outside air can approach the test mass. This gives a 
noise contribution from a single test mass equal to: 



5^(1/1) = 87T^cos{7T[p-l]/2)T{p + 2) 
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(16) 



where x — 27r/rmin/u. 

For typical values v ~ lOm/s and r,„in ~ 5m, at fre- 
quencies above 10 Hz or so, one has x ^ 30 or more, well 
into the exponentially damped regime of the Bessel func- 
tions. Even for gale-force winds of 30m/s or so, the argu- 
ment X of the Bessel functions will still be of order 10 or 
more. The asymptotic expansions oi Kq, Ki, and K2 give 



K^ix) - Ki{x) ^ Kfix) ~ K„ix)Kiix) 



'/{2x'). 



Also, I note that the total noise will be dominated by 
the contribution from the two end stations, which have 
smaller r-n-iin than the corner station. So the total noise 
in the interferometer for uniform airflow is: 



5^(1/1) = (2^)3cos(7r[p-l]/2)r(p+2) 
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I consider "typical" values oi p = 2/3, p = 1.3kg m ^, 
T = 300K, L = 4000m, and r^^^ = 5m. The parame- 
ters V and can vary on a minute-by-minute basis, and 
should really be measured at the site of a given interfer- 
ometer. However, ^ 0.2K^m~^/'^ is a typical daytime 
peak temperature fluctuation index and v ^ 20m/s 
might be typical of a fairly windy day. At frequencies 
around 10 Hz, this gives a noise spectrum of: 



5,(1/1)^ (1.6x10-40Hz-1) 
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lOHz 



-23/3 



10-"(//10Hz) ^ 

(18) 

The two dotted lines in Fig. ^ show the gravity gradient 
noise spectra computed from Eq. ( p^ for wind speeds of 
lOm/s and 30m/s. 

Clearly the exponential cutoff makes this a negligible 
source of noise for LIGO or similar detectors. Physi- 
cally this cutoff arises from the fact that the gravity from 
a particular temperature perturbation passing near the 
end station will affect the test mass coherently over the 
second or so that it takes to travel the width of the end 
station. Thus, even though the temperature noise spec- 
trum has a high-frequency power law tail (reflecting the 
fact that temperature perturbations exist on all length- 
scales), the gravity gradient signal will have this much 
sharper exponentially cut-off tail. 



B. Potential flow near the end station 

As described above, uniform airflow is not likely to 
produce much atmospheric gravity gradient noise in the 
pass-band of interferometric detectors, since the shortest 
timescale over which the gradients change is of order the 
wind crossing time of the interferometer buildings. How- 
ever, if an air pocket could be made to accelerate over 
shorter timescales, it might produce a stronger gravity 
gradient signal at high frequencies. 

One possibility is the acceleration of the air as it is 
forced up and around the wall of an end station: stream- 
lines that approach the ground-wall corner of the end 
station can have curvature scales much shorter than the 
building size. Treating the flow as incompressible and 
vorticity-free, the resulting velocity field near the corner 
is (p. 27 of §): 



Vr, = -~2AX , = 2AZ 



(19) 



where A = v/2R is a constant, and i?, X, and Z are 
measured from the corner. This approximation is clearly 
only good near the corner, since it gives velocity increas- 
ing monotonically with radius; one would expect v to ap- 
proach the free-streaming airspeed y at a distances Rm 
of order the building height. Fig. || shows a schematic of 
the airflow near the corner. 

It is clear from Eq. ( |l9[ ) that, although the streamlines 
are sharply curved near the corner of the end station, 
the advection speed is smaller in direct proportion. The 
shortest timescale over which the motion can change is 
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FIG. 4: Schematic of a uniform airflow streamline passing an interferometer corner station. 9 is the angle between the airflow 
and the axis of the interferometer arm, ro is the distance from the test mass to the streamline at closest approach, ro cos ■)/' is 
the projection of this distance onto the ground, and v =constant is the velocity of the airflow. 




FIG. 5: Schematic of the airflow streamlines around the cor- 
ner of an end station, assuming an incompressible vorticity- 
free flow. The dotted lines are the streamlines, rmin is the 
minimum distance that the airflow can approach the test 
mass, Rm is the scale distance at which the flow velocity v 
approaches the free-streaming speed V , and X, Z, and R are 
coordinates used to describe the velocity field. 



of order Rm/V for all streamlines. Thus the streamlines 
close to the corner will contribute no more high-frequency 
noise than the streamlines further out, at distances of 
order Rm, and the spectrum should not differ greatly 
from the one for uniform flow, Eqs. (p7|), (p^. 



C. Vortices 

Perhaps the most serious contenders for high-frequency 
atmospheric gravity gradient noise are circulating vor- 



tices of air near the end stations. This is somewhat 
stretching the assumptions of the formalism I have estab- 
lished, since I had previously separated the effects of the 
homogeneous turbulence (which establishes the temper- 
ature perturbations on various lengthscales), and wind 
flow (which carries these perturbations past the instru- 
ment). However, the results in this section are only ex- 
pected to be good to order of magnitude anyway, in the 
absence of detailed hydrodynamic analysis of airflow past 
a particular interferometer building. I therefore apply the 
above formalism to a simple model for a turbulent-like 
flow past an interferometer end station. 

For simplicity, I specialize to airflow along the axis of 
an interferometer arm, since these give the largest grav- 
itational accelerations. A simple model for the turbu- 
lent flow is to take a uniform flow and then add a cy- 
cloidal motion to it. This gives x{t) ^ vt — Rsm{vt/R), 
z{t) ^ ro - Rcos{vt/R), r{t) = \/z'^{t) + x'^{t), where 
R is the radius of the cycloidal motion, and ro ~ rmin 
is the distance from the test mass to the unperturbed 
streamline. I treat the speed v along the streamline as a 
constant. If R is also a constant, then one would expect 
the Fourier transform of G{t) = x{t)/r^{t) [Eq. (|l|)] to 
have a spike at frequency / = v/2'kR, with a width of 
order v/tq. However, to give a somewhat more re- 
alistic behaviour for R, I treat it as growing from zero 
at the leading edge of the end station to some scale 
value Rq over the half-length ~ rmin of the end sta- 
tion: R{t) — Roy^vt/raiin- The square-root dependence 
mimics the less-than-linear growth of the thickness of a 
boundary layer. Although quite crude, this model covers 
the essential features of a turbulent flow past the build- 
ing: a uniform translation, accompanied by circulating 
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2R„ 




FIG. 6: Schematic of a vortex-like streamlines across the 
top of an end station. The quantity rmin is the minimum 
distance that the airflow can approach the test mass, and Ro 
is the typical radius of circulation at this distance. 



motions over a range of radii with a scale set by Rq, with 
both the uniform flow rate and the circulation speed set 
by the free-streaming airspeed v. A typical streamline of 
this type is shown in Fig. ^. 

The Fourier transform of G{t) is too complicated to 
perform analytically, but is simple enough to compute 
numerically using a fast Fourier transform. The result- 
ing G(/) has the usual exponential cutoff with frequency 
scale v/2nriaim as for a smooth streamline, but then rises 
to a second peak value of ~ (4/i')\ARo/^^iiii ^ 
quency / ~ 0.06w/i?o and decreases from there as /"^. 
This high-frequency tail is the power law that one would 
expect from the cusps on the bottom of the cycloid. Since 
I treat v as constant, F{f) = vP~^^G{f). The cross- 
sectional area of streamlines that contribute significantly 
around this peak frequency is of order ~ 2rmin-Ro- Plug- 
ging these into Eq. (10) one gets a noise spectral density 
of 




FIG. 7: Plot of dimensionless noise amplitude yfShiJ) 
versus frequency /. The solid line is the projected noise floor 
for advanced LIGO detectors. The two dotted lines are the 
gravity gradient noise levels caused by temperature perturba- 
tions advected along smooth streamlines at lOm/s (left curve) 
and 30m/s (right curve). The dashed lines are the noise levels 
caused by temperature pertur batio ns advected along cycloidal 
vortices, as described in Sec. Ill C| , at lOm/s (left curve) and 
30m/s (right curve); in each case the size scales Rq of the vor- 
tices have been tuned to maximize the noise at 10 Hz. Even 
given these fine-tunings, the gravity gradient noise curve in 
the worst-case scenario is still nearly an order of magnitude 
(in amplitude) below the advanced LIGO sensitivity. 



(max) 



1.5 X 10^ 



Gp_ 
LT 



max mm 



(20) 



at the peak frequency of / ~ 0.06ti/i?0j where I have 
assumed p = 2/3. 

For wind speeds around lOm/s, a cycloid radius i?o ^ 
0.06m(z)/10m/s) puts the noise peak at the 10 Hz seis- 
mic noise wall for advanced LIGO detectors. The atmo- 
spheric gravitational noise contribution from a single end 
station is then: 



?Mmax) ~ (1.3 X 10-49HZ- 
20/3 



0.2Km 



-2/3 



lOm/s 



5m 



(21) 



This is over five orders of magnitude below the expected 
advanced LIGO noise floor of 2 x lO^'^^Hz"^ at 10 Hz. 
Gale-force winds (w ^ 30m/s) will bring this up to 
2 X 10~^^Hz~^, still an order of magnitude below the 
advanced LIGO noise curve. The dashed lines in Fig. |^ 
show the actual data from the numerical Fourier trans- 
forms for these two cases. 



Since even the worst-case estimate is still an order of 
magnitude below the advanced LIGO noise floor (as well 
as the seismic gravity gradient noise floor in |3|), it seems 
unlikely that turbulent vortices will be sufficient to raise 
the atmospheric gravitational noise to significant levels, 
even given the approximations made in this analysis. 
However, to settle this matter definitively would require 
a much more sophisticated numerical analysis of the tem- 
perature perturbations and airflow past the buildings of 
a particular facility. 



IV. SHOCKWAVES 

Although atmospheric pressure waves are unlikely to 
be a significant source of gravity gradient noise in inter- 
ferometric gravitational- wave detectors, the sudden pres- 
sure changes caused by atmospheric Shockwaves could 
potentially produce detectable transient signals in the 
detector, if such shocks occur in the vicinity of the de- 
tectors. Shocks are specifically a matter of concern be- 
cause they can produce significant pressure changes over 
timescales less than 0.1 s, corresponding to the lower end 
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of the pass-bands of most interferometric detectors. 

Consider a shock that produces a sudden jump in air 
pressure in the vicinity of one of the interferometer test 
masses: Sp{t) — ApQ{t — to). It is a simple matter to 
take the Fourier transform and apply Eq. to obtain: 

Hf) = o fTf4 — cos(0)C(2^/r„,i„/c) . (22) 

Here, 9 is the angle between the interferometer arm and 
the normal to the shock front. If the shock has a finite rise 
time At, we can mimic this analytically by convolving in 
the time domain with a Gaussian of width a ^ At. In 
the frequency domain this multiplies our amplitude by 
a Gaussian of width a ^/At, giving an exponential 
cutoff at frequencies above 1/At. Typical shocks from, 
for instance, supersonic objects have rise times on the 
order of a few milliseconds corresponding to cutoff 
frequencies of a few hundred Hz. However, one expects 
the dominant contribution of the signal to come from 
much lower frequencies, before the building-size cutoff 
factor C kicks in. 

Shocks are transient phenomena that will produce sig- 
nals in the detector, rather than raising the noise floor. 
What one would like to know is what signal-to-noise ratio 
the shock will produce. In general this depends on what 
filters one is using to search for signals, and how well 
these filters overlap with the signal produced by a shock. 
However, the signal from a shock is likely to overlap quite 
well with templates designed to search for generic impul- 
sive phenomena; such templates are likely to be used in 
advanced interferometers as control over nonstationary 
instrumental noise improves. Thus it is reasonable to 
consider the signal-to-noise ratio p that a matched filter 
would give |9|: 



Mf)h*if) 
Shif) 



df 



\vH.f)\' 

fSnif) 



din/. (23) 



This shows, roughly speaking, that the relative magni- 
tude of the dimensionless signal amplitude 2\fh{f)\ and 
the dimensionless noise amplitude \/ f Sh(f) over a loga- 
rithmic frequency interval gives a good indication of the 
signal-to-noise ratio produced in the detector. 



A. Sonic booms 

Sonic booms caused by supersonic bodies are one ex- 
ample of atmospheric shocks that might affect interfer- 
ometric gravitational-wave detectors. Direct Shockwaves 
from a supersonic aircraft will typically hit the ground 
in a "carpet" about 15-20 km wide under the aircraft's 
flight path. Outside this carpet, the temperature gradi- 
ent near the ground will completely reflect the Shockwave 
before it touches down. However, Shockwaves will also 
reflect downward off of the temperature inversion in the 
stratosphere and thermosphere, forming secondary and 



higher-order "carpets" out to many hundreds of kilome- 
tres [|o| . The presence or absence of these higher-order 
waves can depend quite sensitively on conditions in the 
upper atmosphere. 

A detailed prediction of these effects is beyond the 
scope of this paper. However, to give an indication of 
their potential seriousness, I will consider what would 
happen if a supersonic aircraft were actually to overfly 
the instrument at a height of several kilometres. 

The characteristic profile of a sonic boom is a sym- 
metric N-wave, consisting of a shock that increases the 
pressure by an amount Ap, followed by a smooth de- 
crease in pressure of 2Ap over a time At, followed by a 
second rising shock Ap to restore the ambient pressure. 
According to Eq. (9.78) of ||ll|, the strength of the shocks 
is: 

— ^ I '^T^ln iM' - , (24) 

p (7-fl)^/^ 

where 7 is the adiabatic coefficient of air (1.4 at normal 
temperatures) , M is the Mach number of the aircraft (its 
speed divided by the sound speed), k is a dimensionless 
form factor that depends on the shape of the aircraft 
(typically around 1), I is the length of the aircraft, 8 is 
the ratio of the aircraft's typical thickness to its length, 
and r is the closest distance that the aircraft came to the 
point of measurement. Between the two shocks, the rate 
of change of pressure in a direction e^, parallel to the line 
of flight is given, in Eq. (9.80) of (ill, as: 



\dp ^ 7 (M2 _ 1)1/2 
p dx 7 -|- 1 r 



(25) 



The shock fronts move outward at the sound speed c in 
the direction orthogonal to their surface, while the entire 
cone travels in the direction along with the aircraft 
at a speed Mc. The total change of pressure between 
the two shocks is 2Ap. From these facts and Eqs. ( |2^ ) 
and (p5|), one can show that the time between the two 
shocks is: 



At « 2^/4(7 + 1)1/2 



M 



(M2 - 1)3/8 



;3/4 1/4 

kS . (26) 



On frequency scales higher than 1/At (typically a few 
Hz for a supersonic aircraft a few kilometres away), the 
sonic boom looks like simple Heaviside shocks, giving 
Sp{f) ^ 1/f. At lower frequencies, though, the entire N- 
wave looks like the derivative of a (5-distribution, giving 
^Pif) ~ /• Performing the Fourier transform analytically 
and plugging into Eq. (0), one obtains: 



kf) = T^i— cos(0)C(27r/r,„i„/c) 



Att^jL p 
"sin(7r/At) 



TTfAt 



cos(7r/At) 



(27) 



where to is the time when the midpoint of the N-wave 
crosses the detector. As expected, the amplitude goes 
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roughly as /~^, except for frequencies less than 1/At, 
where it goes as f~^. 

Now let us plug in some typical numbers. The numbers 
G = 6.67 X IQ-iim^kg^^s-^ p = LSkgrn^^, 7 = 1.4, 
c ~ 332m s^^, and L ~ 4000m can be treated as con- 
stant. A supersonic jet aircraft might have a length of 
I = 10m, a typical diameter of dl — 2m, and be traveling 
at Mach M = 1.5 at a distance of r = 10km or so. Let 
cosO be 1 for an upper limit. Then At ^ 0.2s, and for 
frequencies / > lOHz the dimcnsionless signal amplitude 
is: 



2|W)I 



1.4 X 10-i9(Af2 _ l)i/8c(2^/r„,i„/c) 



10m 



\ 3/4 

) ( 10km) 



-3/4 



lOHz 



_S_ 

02 



(28) 



This is three orders of magnitude above the expected 
noise floor of ^/J%JJ) - 1.4 x lO'^^ at 10 Hz for ad- 
vanced LIGO interferometers! 

By contrast, consider a .30-calibre rifle bullet (/ « 
0.025m, 6 w 0.3) passing at Mach 3 within 10m of an 
interferometer test mass. (This stretches the assumption 
of a plane-wave shock front at the test mass, but the or- 
der of magnitude should be correct.) The bullet produces 
a much stronger double shock, but with a time interval 
At ~ 0.5ms, so l/At^ 2kHz. The low-frequency tail of 
this signal will have dimcnsionless amplitude: 



, ~ 23M2(Af2_l)-5/8 

2\fhif)\ ^ 1.8 X 10-23—^ ^ 



0.3 



I 



0.025m 



2.5 

9/4 



/ r N -1/4 / / 



(lOm) 



lOHz 



C(27r/r,„i„/c) 



(29) 



This is nearly an order of magnitude below the dimcn- 
sionless noise amplitude in advanced LIGO, and therefore 
too small to be of any serious concern. 

Fig. H shows more complete gravity gradient signal 
spectra computed using Eq. ( p7| ) with the above param- 
eters for a supersonic aircraft and rifle bullet. These are 
plotted along with the anticipated dimcnsionless noise 
amplitude for advanced LIGO detectors. 
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FIG. 8: Plot comparing dimensionless signal strengths 
2\fh{f)\ of gravity gradients from sonic booms with the di- 
mensionless interferometer noise amplitude yfShiJ), as a 
function of frequency. The solid line is the projected noise 
amplitude for advanced LIGO detectors. The dotted line is 
the signal from a .30 calibre bullet traveling at Mach 3 within 
10m of the end station; the dashed line is the signal from 
a 10m long aircraft traveling at Mach 1.5 at a distance of 
10km. The aircraft's sonic boom can produce a detectable 
gravity gradient signal. 



Alternatively, if the same shock proflle is detected in 
an array of at least three sound sensors, then one can 
determine the direction of propagation of the shock and 
predict the actual induced test-mass motions. The spu- 
rious signal could then be subtracted out of the data 
stream. This is a much trickier procedure, and would 
only be necessary in the unlikely event that signiflcant 
amounts of data were being corrupted. 

In either case, it is clear that infrasound sensors will be 
important environmental monitors for advanced interfer- 
ometric detectors. 



B. Vetoing Shockwave signals 



HIGH-SPEED OBJECTS 



While atmospheric Shockwaves are a potential source 
of spurious signals in gravitational-wave detectors, they 
are easy to veto using environmental sensors. One need 
simply place infrasound microphones outside the build- 
ings and test-mass vacuum enclosures. If these sensors 
detect a pressure change of more than a millibar over 
timescales of 50-100 milliseconds, then one might expect 
spurious signals with dimensionless amplitude of 10~^^ 
in the 10-20 Hz frequency range. The stretch of data 
containing the potential spurion can then be discarded. 



Another potential source of spurious signals in the in- 
terferometer is the gravity gradient caused by the motion 
of an individual massive object past the interferometer 
end station, or the collision of such an object with the 
end station. The latter is particularly serious, since the 
sudden deceleration of the object can produce a signal 
at high frequencies. The issue of human-generated grav- 
ity gradient noise has been addressed in ||^, but there 
are other sources outside the facility that must be con- 
sidered, such as stray bullets and wind-borne debris. In 
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particular, the Hanford LIGO facility is plagued by tum- 
bleweeds, which can produce non-negligible gravity gra- 
dient signals. 

The general formula for the spurious gravitational- 
wave signal produced by a moving object is: 



kf) 



GM 



2-^/* dt , 



(30) 



where M is the mass of the object, r{t) is its distance 
from the test mass as a function of time, and x{t) is its 
distance from the end mass in the direction parallel to the 
interferometer arm. For an object traveling parallel to 
the ground in a straight line at speed v, Eq. (|3^) becomes 
[as in Eq. (Rl)]: 



Hf) = 



GM 
Lv^nf 

-iKi 



27r/r„ 



cos6' 



27r/r 



sin 6 cos ip 



,(31) 



where 9 is the angle between the line of motion and the 



interferometer arm, 



is the distance of closest ap- 



proach between the object and the test mass, V is the 
angle projecting this distance onto the ground, and 
is the time of closest approach. Kq and Ki are modi- 
fied Bessel functions of the second kind of order and 1, 
respectively. 

Under moderately windy conditions (wind speeds up 
to 15m/s or so), tumbleweeds at the Hanford LIGO 
facility will bounce along the ground at 5-lOm/s. In 
stronger winds, the tumbleweeds become airborne, with 
speeds approaching the wind speed; they can fly right 
over the LIGO buildings, or impact with considerable 
force. The same may be true of wind-borne debris at 
other interferometer facilities. However, the value of 
fmi-a is usually at least 5m, so for frequencies above 
10 Hz one has 27r/rinin/i' ^ 10 even for very strong 
winds [v ^ 30m/s). In this regime the Bessel func- 
tions are exponentially damped, so these objects will 
not produce significant gravity gradient signals simply 
by blowing past the instrument. A rifle bullet, on the 
other hand, might be moving around lOOOm/s, putting 
us in the small- argument regime of the Bessel functions, 
where Kf){x) ^ ~\ti{x) and Ki{x) ^ . Taking the 
most dangerous geometry 9 — tt/2, ip — 0, and assuming 
a bullet mass of around 5 grammes, this gives a dimen- 
sionless signal amplitude near 10 Hz of: 



2|A(/)I 



1.6 X 10 
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5g/ \ V 
/l0Hz\ 
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(32) 



This gives a signal-to-noise ratio of about 1 at 10 Hz, a 
bit below the detectable threshold. In fact, even if one 
fine-tunes the bullet speed v, the largest signal-to-noise 
ratio that one can get at 10 Hz is about 2, for a speed of 



around 250m/s. Since events with signals less than about 
5 times the noise will probably be ignored in any case, 
one can conclude that objects flying past a test mass are 
not likely to be serious sources spurious events. 

If an object does not pass smoothly by the interferom- 
eter but instead collides with an end station, the signal 
at ~ 10 Hz can be large even for slow-moving objects: 
it is the deceleration time, not the end-station-crossing 
time, that sets the frequency scale of the signal. Suppose 
an object collides end-on with the end station at a speed 
V, coming to a stop within a distance d at constant accel- 
eration. Let t — denote the time that the object comes 
to rest. The motion of the object is then given by: 



r{t) = xit) 



d ~ vt 
{vtf/Ad 



t < -~2d/v 
-2d/v <t <0 
t > 



(33) 

The Fourier transform oi x{t) / [t] — l/r^(t) is tricky to 
do analytically, so I have relied on numerical fast Fourier 
transforms, and then made approximate analytic fits to 
the result. However, one can qualitatively predict the 
shape of the signal in frequency space. The function 
l/r^(t) starts out near zero and then slowly rises over 
a timescale rmin/i' to a value rf^^y^, then quickly lev- 
els off at that value over a timescale d/v. So on fre- 
quency scales <C f/rmin, l/'^^(^) looks like a step func- 
tion, whose Fourier transform goes as 1//. On frequency 
scales > w/rniin but <C f/rf, one sees the deceleration as 
a cusp (discontinuous first derivative), giving a Fourier 
transform that goes as 1//^. On frequency scales > v/d, 
the deceleration appears smooth, but the onset of deceler- 
ation in sudden, giving a 1/ behaviour. The following 
gives a good fit to the numerical Fourier transform: 



oo g27ri/t 
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■ dt 
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(34) 



The two breakpoints separating the three branches are 
at frequencies /i — OAv/r^nin and /2 = 0.24w/d. More 
precisely, since the deceleration occurs over a well-defined 
time 2d/v, the third branch of the Fourier transform is 
oscillatory with nodes every v/2d in frequency space; the 
functional fit in Eq. ( |3^ ) is an envelope containing these 
oscillations. 

For a 5g bullet striking an end station at lOOOm/s, 
the signal at 10 Hz is dominated by the low-frequency 
tail regardless of stopping distance. The dimensionless 
signal amplitude is then: 



2|/M/)|'-3xlO-^^ — 
5g 



M\ /lOHz\ / 5m 



(35) 



The signal amplitude is about the same as for a pass- 
ing bullet, although the dependence on / (and r^i^ and 
v) is different. For a tumbleweed or other wind-borne 
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object, by contrast, the sudden deceleration can create 
significant high-frequency noise. A typical tumbleweed 
at Hanford has a mass of AI = 0.1kg and a diameter of 
0.4m, and can be compressed by about half that amount 
{d — 0.2m). Larger weeds can be twice as large in di- 
ameter, putting their masses in the 0.5-lkg range [^ . 
For moderate to high wind speeds {v — 10-30m/s), the 
signal at frequencies above 10 Hz is dominated by the 
second branch of Eq. (^) , giving a dimensionless signal 
amplitude of: 



2|/M/)I - 5 X 10-21 



M 
1kg 



/10Hz\' 



lOm/s 

3 



5m 



(36) 



Thus even a "typical" 0.1kg weed at these speeds will 
produce a signal-to-noise ratio of around 4 at 10 Hz in 
advanced LIGO interferometers, which is in danger of 
being interpreted as a real gravitational- wave event. In 
a more extreme case, a 1kg tumbleweed blown airborne 
by a strong 30m/s wind will produce a signal 100 times 
higher than the noise at 10 Hz, which is easily detectable. 

Fig. ^ shows the signal spectra for the objects dis- 
cussed above, plotted against the dimensionless noise 
amplitude expected in advanced LIGO detectors. Since 
tumbleweeds are a potential source of spurious detectable 
events, one should consider ways to reduce the tumble- 
weed gravity gradient noise. Fortunately, the signal goes 
as ?'~fn, so a simple fence preventing the weeds from ap- 
proaching the end station should be sufficient. A fence 
30m out from the end station will reduce the signal-to- 
noise ratio to 1 for tumbleweed masses up to 1kg and 
speeds up to 30m/s, reducing the risk of spurious events. 



VI. CONCLUSIONS 




FIG. 9: Plot comparing dimensionless signal strengths 
2\fh{f)\ of gravity gradients from airborne objects with the 

dimensionless interferometer noise amplitude \/ fSh{f), as a 
function of frequency. The solid line is the projected noise 
amplitude for advanced LIGO detectors. The short-dashed 
line is the signal from a .30 calibre bullet passing next to 
the interferometer end station at lOOOm/s, while the dotted 
line is the signal from that same bullet colliding with the end 
station. The long-dashed line is the signal from a large 1kg 
tumbleweed passing next to the end station at 30m/s, while 
the dot-dashed line is the signal from that same tumbleweed 
colliding with the end station, assuming that it compresses 
by 20cm on impact. The tumbleweed mass and speed are 
both near the upper end of the expected range, representing 
the greatest danger of spurious signals in the detector; evi- 
dently these tumbleweed signals can be detected quite easily 
by advanced LIGO detectors. 



This paper has studied two sources of background grav- 
ity gradient noise, from infrasonic atmospheric pressure 
waves and from wind-advected temperature perturba- 
tions, in order to determine whether they constitute a 
limiting noise floor for interferometric gravitational-wave 
detectors — in particular, for the "advanced" LIGO detec- 
tors projected in . The paper also analyzed two sources 
of gravity gradient signals, from transient atmospheric 
Shockwaves and from massive airborne bodies, to deter- 
mine whether they would constitute detectable spurious 
events in these interferometers. The following summa- 
rizes the results and suggests possible further work that 
may need to be done. 

Current estimates suggest that infrasonic pressure 
waves will not be a significant source of gravity gradi- 
ent noise, being over two orders of magnitude below the 
advanced LIGO noise floor at 10 Hz. Nonetheless, these 
estimates are not based on actual noise measurements 
at an interferometer site, so infrasound measurements at 
these sites are recommended to confirm them. Further 



empirical studies might also analyze the specific effects 
of building shapes and of infrasound coherence lengths 
on the noise spectrum, although these refinements would 
likely only serve to reduce noise estimates above 15 Hz 
or so. 

Wind-advected temperature perturbations, although 
the dominant source of atmospheric density fluctuations, 
do not produce significant high-frequency gravity gradi- 
ent noise, due to the long times that any particular pocket 
of warm or cool air spends in the vicinity of an interferom- 
eter test mass. A possible exception is when the airfiow 
forms vortices around the interferometer buildings, since 
this will produce a noise spectrum peaked around the 
typical vortex circulation frequencies near the test mass. 
The current crude analysis of these effects suggests that 
the noise is still an order of magnitude below the ad- 
vanced LIGO noise floor at 10 Hz even in the worst-case 
scenarios, but the model could be improved significantly. 
Numerical models of the airfiow and of temperature per- 
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turbations near an interferometer building may be re- 
quired to settle this issue definitively. 

Gravity gradients from atmospheric Shockwaves are 
potentially serious sources of spurious signals in interfer- 
ometric gravitational-wave detectors. For instance, the 
sonic boom from a supersonic aircraft overflying an ad- 
vanced LIGO detector could produce signal-to-noise ra- 
tios of several hundred. Although such overflights are 
expected to be rare or nonexistent, they point out the po- 
tential seriousness of shocks from weaker or more distant 
sources, even if the signals are several orders of magni- 
tude smaller. It is therefore strongly recommended that 
advanced interferometric detectors include infrasonic de- 
tectors as environmental monitors. Such sensors could 
easily be employed to veto spurious atmospheric gravity- 
gradient events. 

Gravity gradients from wind-borne objects such as 
tumbleweeds are another possible source of spurious 
events in gravitational-wave detectors, if these objects 
are allowed to collide with the buildings housing the in- 
terferometers. Fences or other structures should be used 
to keep these objects at least 30 metres from the test 
masses, in order to eliminate the risk of spurious signals. 

Obviously the number of things that can affect inter- 
ferometric detectors through gravity gradient forces is 
immense; I have considered here only the few sources 
that I considered the most worrisome. I encourage other 
researchers to consider the implications of this often- 
neglected effect. 
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APPENDIX A: THE TEMPERATURE NOISE 
SPECTRUM 

This appendix presents a more rigorous mathematical 
derivation of Eq. ( |l0| ) used in Sec. III. 

Consider a time-varying field of temperature pertur- 
bations 6T{r, t) about some average temperature T. 
This produces a gravitational perturbation gxit) — 
J dV Gpxr~^{ST/T), where p is the average air density, 
and the x-axis is along the interferometer arm. The spec- 
tral density of gravity gradient noise 5*9(1/1) is given by 
twice the Fourier transform of the gravity autocorrelation 
Cg{T) = (g{t)git + T)). Thus: 



{ST{r,t)ST{r\t + T)}e 



2771 fr 



(Al) 



The temperature noise measured at a point Tq, on the 
other hand, is given simply by: 



sAlfl) 



dT(ST{ro,t)ST{ro,t + T))e 



27TifT 



(A2) 



On sufficiently small scales, the temperature pertur- 
bations in the Earth's turbulent boundary layer can be 
treated as homogeneous and isotropic. The expected 
squared temperature difference between two points is 
then a function only of their separation: {[ST{r) — 5T{r + 
Ar)]2) = DrdlArll). The function DriAr) is called the 
temperature structure function of the atmosphere, and 
for small Ar reduces to a power law Drir) = c^Ar^. 
If a wind with speed v blows these perturbations past 
a measuring station, the temperature autocorrelation is 
{ST{ro,t)ST{ro,t + T)) ^ - (l/2)4(wr)P for smaU r, 
where is the mean squared temperature fluctuation. 
This results in a high-frequency power law tail: 

Srilfl) = 4t;^'(27r/)-(^'+i)r(p + 1) sin(p^/2) . (A3) 

Turbulent mixing theory, as well as micrometeorologi- 
cal measurements of 5t(|/|), show that the value of p 
is normally 2/3, characteristic of a type of turbulence 
known as Kolmogorov turbulence. See, for example, ||] 
for discussion of this type of turbulent mixing, also ||] 
and references therein. 

We are interested in Sg{\f\), which is somewhat trick- 
ier to calculate than '5't(|/|), since it involves a correla- 
tion between points separated in space as well as time. 
However, chaotic turbulence will almost certainly de- 
stroy high-frequency correlations between widely sepa- 
rated points, so the high-frequency behaviour of Sg{\f\) 
will come from correlations between nearby points. That 
is, the high-frequency support of {5T{r,t)ST{r' ,t + t)) 
will come from those points r' at time t + t whose fluid 
elements were near r at time t. 

Consider two fluid elements moving along paths S and 
S' passing through r and r', respectively. This is shown 
schematically in Fig. 11^. Let ro be the distance from r to 
the nearest point on S^^, and let tq be the time it would 
take a pocket of air at r' to be carried to this point on 
S'. In order for these points r and r' to contribute to 
the high-frequency component of the spectrum, the dis- 
tance ro must be fairly small, of the order v/f where 
V is the wind speed past r. I treat the streamlines as 
relatively straight over these scales, in which case the 
temperature perturbations at (r, i) and (r',f-t-r) corre- 
spond to physical pockets of air separated by a distance 
\/rQ + u(r)2(r — tq)^, for r near tq. Assuming that the 
T-dependence of the correlation function is due entirely 
to this advection, the correlation function can be written 
explicitly as: 



SA\f\) = 



dr dV dV' 



\r'Y 



{ST{v, t)6T{v', t + T))=al-^ [rl + v\t ■ 



To 



(A4) 
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where typically p k, 2/3. 

This term contains the entire r-dependence of the grav- 
ity perturbation in Eq. (Al), so the first integral I do is 
the Fourier integral over r. This integral has the form 
!T^{fi'^+x^Y-^/'^e"'''' dx = 2i:-^/'^{2(3/aY cos{-Kv)T{v + 
l/2)K_^{a(3). Formally the integral diverges for > 1/2, 
but the closed-form expression remains approximately 
correct for large a provided the integrand is cut off 
smoothly for large x 3> 1/a. Physically this corre- 
sponds to the fact that a smooth, larg e-scale cutoff of the 
temperature correlations in Eq. (A4) will not affect the 
high-frequency component of the temperature noise. For 
horizontal winds near the ground the spatial correlation 
function is cut off on horizontal distance scales of ^ 50 
times the fluid elements' altitude z above ground, giv- 
ing a low-frequency cutoff around ~ 0.02u/z (Fig. 1.A4 
of [^). Typically this will be below the relevant fre- 
quency range for intcrferometric detectors; I ignore it 
to obtain pessimistic (upper-limit) noise estimates. The 
high-frequency noise tail is then: 



Sm) = (^) JdvJ dV 

X 4z;P(V27r/)-(f+i)a(p) 
X i^(p+i)/2(27r/ro/t;)e2^^^^° 



^3 ^ j^^ ) ^ 



(A5) 



where a(p) = -27r-i/2 cos(7r[p+ l]/2)r([p-h2]/2) w 0.873 
for p = 2/3. 

Next is the integral over dV' . The exponential decay 
of the Bessel function if(p+i)/2 restricts the support of 
this integral to values rg ^ t;/27r/, representing the fact 
that high-frequency fluctuations can only arise from the 
rapid change of the spatial correlation function over small 
lengthscales. This range in tq defines a narrow bundle of 
streamlines S' about the streamline 5* passing through 
r. I assume that the size of this bundle is less than the 
distance from the test mass to the bundle, so that the 
values of x' and r' on a given S' can be replaced with the 
nearby values on S. Now for points near r the volume el- 
ement dV can be written in terms of the new parameters 
To and To as dV — 2TTrQ dro v(r)dTQ. Since the airflow, 
being very subsonic, is nearly incompressible, the volume 
element retains this form for all points along the bundle: 
if v{r') decreases below w(r), for instance, the length el- 
ement v{r')dTQ will decrease, but the cross-section of the 
bundle (i.e., the relevant range of 2TTrodro) will increase 
to compensate. Plugging in this volume element and in- 
tegrating over ro, one obtains: 



XX 



2 



^4«^(7r/)-(P+i) 
(A6) 



27ri/To 



X a(p)2^r(b + 3]/2) [ — I e 



Let t' be a new time coordinate denoting the time it 
takes for an air pocket to reach r' from some fixed refer- 
ence plane that crosses all streamlines orthogonally, and 
t be the corresponding coordinate for the point r. Then 
To is just t' — t, and the volume element dV can be written 
as w dt dA, where dA is a cross-sectional area element on 
the reference plane, and w is the wind speed across that 
area element. Plugging this in, and ignoring the spatial 
separation between S and S", one obtains: 

= (^)'4(^/)-(^+=^'«(p)(7r/2)r(b + 3]/2) 



X / wdA 

l{S} 



±vP+h-^^'f'dt 



x' 



(A7) 



where I have used J^^y to denote an integral over the 
entire reference plane; i.e., over all streamlines S. 



15 



To 



s ^■ 



v(r) 





{S} 



FIG. 10: Schematic showing the parameters used to describe streamlines of air flowing past an interferometer test mass. {S} 
denotes a plane intersecting all streamlines; S and S' are two such streamlines passing through the points r and r', respectively, 
where the coordinate origin is centered on the test mass k. The distance between the streamlines, measured at r, is ro. The 
positions of the points r and r' along the streamlines S and S' are parameterized by the times t and t' that it would take for 
a pocket of air to move from {S} to r or r'; to is the difference t — t' . v(r) is the wind velocity along S through the point r; w 
is the wind velocity through the plane {S}. 



The noise in the gravitational-wave signal h{t) due 
to the gradients at a given test mass is Sh{\f\) — 
{27rf)^'^Sg{\f\)/L'^, where L is the interferometer arm 



length. The noise at each test mass adds incoherently to 
the total signal. Combining these with Eq. (A7) yields 
the result given in Eq. ([lO|). 
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